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We present a Raman spectroscopic investigation of the Herbertsmithite ZnCu3(OH)6Cl2, the first 
realization of a Heisenberg s = 1/2 antiferromagnet on a perfect kagome lattice. The magnetic 
excitation spectrum of this compound is dominated by two components, a high temperature quasi 
elastic signal and a low temperature, broad maximum. The latter has a linear low energy slope and 
extends to high energy. We have investigated the temperature dependence and symmetry properties 
of both signals. Our data agree with previous calculations and point to a spin liquid ground state. 
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I. INTRODUCTION 

Strongly frustrated spin systems are in the focus of 
many recent works as the competition of exchange in- 
teractions leads to interesting ground states and a rich 
low energy excitation spectrum. Of these systems, the 
s = 1/2 Heisenberg antiferromagnet on the kagome lat- 
tice has been of special interest, since it displays a unique 
quantum disordered and highly degenerate ground state, 
i.e. a spin liquid^. A kagome lattice consists of a 
two-dimensional (2D) arrangement of corner-sharing tri- 
angles. In systems with large spin so called weather- 
vane modes, low energy collective rotations of spins on 
a hexagon, destabilize any order. For s = 1/2 a quan- 
tum superposition of an infinitely dense spectrum of low 
energy singlet and triplet excitations are supposed to ex- 
ist. Defects on the kagome lattice may lead to a gapless 
valence bond glass phased with enhanced singlet bonds 
opposite to the defects. Together with antisymmetric 
Dzyaloshinskii-Moriya (DM) interactions a quantum crit- 
ical point is induced^. Despite these and other— theoreti- 
cal investigations, the understanding of the ground state 
properties and the excitation spectrum of the quantum 
kagome system is not completely satisfactory. 

Earlier experimental investigations of compounds that 
are close in realizing the kagome layer are hampered ei- 
ther by interlayer interactions, a non-uniform, more com- 
plex distribution of exchange coupling constants or spin 
anisotropics. The first reported material with 2D kagome 
layers was SrCr8Ga40i9^ with s = 3/2. It shows a spin- 
glass transition for temperatures below T = 3.3 K despite 
the much larger interaction strength. In the s = 3/2 com- 
pound Yo.5Cao.5BaCo407^ short ranged chiral spin cor- 
relations develop at low temperature, while the s = 1/2 
compound RbCusSnFiT^ experiences a spin gap that sep- 
arates the singlet groundstate from low energy triplet ex- 
citations. 



Presently the compound ZnCu3(OII)6Cl2 is the only 
well characterized realization of the s = 1/2 Heisen- 
berg antiferromagnet on a perfect kagome lattice^i^. It 
exists in two structurally slightly different forms: a) 
Kapellasitei^ and b) Herbertsmithite. Small single crys- 
tals, however, are up to now only available as Herbert- 
smithite. In this compound, the Cu^+ ions form the tri- 
angular kagome lattice in the crystallographic ab-layer 
and are coupled by oxygen atoms via superexchange 
pathways, leading to strong antiferromagnetic correla- 
tions of about J w 180 K. Despite this strong coupling, 
susceptibility measurements reveal no magnetic ordering 
down to 50 mKii. 

The planes of these kagome structures are separated by 
intermediate layers of Zn and CI atoms and can therefore 
ideally be considered as 2 dimensional layers with neg- 
ligible magnetic interactions along the third dimension. 
However, Zn/Cu antisite disorder of at about 6%i^ is be- 
lieved to be present in previously studied Herbertsmithite 
samples. Hence, the presence of a small concentration of 
impurities may affect the spin correlations in the kagome 
plane. 

Muon spin rotation (^SR) measurements on the parat- 
acamite Zno.5Cu3.5(OH)gCl2 showed a partly frozen 
ground statoH, which however disappears for higher 
Zn/Cu ratio. For the Herbertsmithite, i.e. a Zn to Cu 
ratio of 1/3, /iSR did not reveal any ordering down to 50 
mK. Electron spin resonance measurements revealed DM 
interactions of about D = 0.08 • J magnitud&ii, which 
causes a ferromagnetic-like increase of susceptibility at 
low temperature. These and other experiments have been 
performed either on natural (mineral) crystals or powder 
samples^. Only recently, Herbertsmithite single crystals 
were successfully grown via hydrothermal synthesis^. 

Our approach of investigating the excitations in the 
kagome lattice is using inelastic light scattering (Raman 
scattering). This is a unique method in the sense that it 



is sensitive to low energy topological, electronic, or mag- 
netic excitations. Singlet modes (s = 0) can be observed, 
which may otherwise only be probed in specific heat mea- 
surements. In particular, it has recently been suggestecU^ 
that studying the polarization dependence of the mag- 
netic Raman scattering contribution should allow to dis- 
tinguish between different possible ground states such as 
a valence bond crystal or a spin liquid. Of further inter- 
est is the spectral distribution of magnetic scattering as 
characteristic low energy slopes or sharp modes may be 
observed as signs of the topological character of the spin 
correlationsiii^. 



II. EXPERIMENTAL 

This study deals with two different types of samples: 
(i) Natural crystallites from mineral sources, which show 
a dark green color and are slightly transparent. A se- 
lected crystallite from this source with the approximate 
volume of about 0.5 mm'^ was studied via Buerger pre- 
cession methodic to confirm the sample's single domain 
structure and to uncover a natural surface that is parallel 
to the crystallographic ab-plane, in which the kagome lat- 
tice is realized. Its stoichiometry, Zno.8Cu3.2(OII)6Cl2, 
is slightly shifted from optimum, (ii) A hydrothermally 
synthesized single crystal with a size of about 1 x 0.2 x 
0.2 mm'^, which shows a very regular habitus, was stud- 
ied in addition. It appears in light green color and high 
transparency and has a stoichiometry very close to the 
perfect kagome, i.e. ZnCu3(OII)6Cl2. 

Raman spectroscopic studies were performed in quasi- 
backscattering geometry with a solid state laser (A = 532 
nm) and 1 mW laser power. Room temperature measure- 
ments were carried out under ambient conditions, while 
low-temperature measurements were done in an evacu- 
ated, closed-cycle cryostat. Investigations of the polar- 
ization dependence were performed by rotating the sam- 
ple within the ab plane. In the 0° orientation, the a axis 
was approximately perpendicular to both the polariza- 
tion and the k vector of the incident laser light. The 
spectra were collected via a triple spectrometer (Dilor- 
XY-500) by a liquid nitrogen cooled CCD (HORIBA 
Jobin Yvon, Spectrum One CCD-3000V). In a previous 
experiment powder samples of Zn2;Cu4_2;(OH)6Cl2 (with 
X = 0.5, 0.85, 1.0, 1.1, 1.2 and 1.4) were investigated to es- 
timate the stoichiometry of the single crystals by compar- 
ing their phonon frequencies^^. These samples, however, 
do not allow to successfully study the weaker magnetic 
or topological Raman scattering in the Herbertsmithite. 



III. RESULTS 

Figure 1 a) shows Raman spectra of the natural and 
of the synthesized single crystal, measured at room tem- 
perature (T = 295 K) and in xx polarization. Rotating 
the samples within the crystallographic ab plane (with 
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FIG. 1. a) Raman spectrum of natural and synthesized single 
crystal at room temperature in xx polarization, b) left: inten- 
sity of the quasi elastic signal at RT as function of the angle 
between the light scattering polarization and a fixed direction 
within the ab plane as described in the text. The rotational 
anisotropy of this scattering intensity is consistent with that 
of Aig phonon modes; right: intensity of the Eg mode at 363 
cm^^ plotted in a similar fashion, c) Comparison of phonon 
modes in ZnCu3(OF)6Cl2 assigned to H/D vibrations com- 
paring a hydrated (y=H) and a deuterated {¥=!)) powder 
sample. 



the axis of rotation along the k vector of the incident 
light) allows to distinguish between modes of Aig and Eg 
symmetry representation. Furthermore, the previously 
discussed rotational anisotropy of the magnetic modes 
can be probed. 

We observe 7 sharp modes at low to intermediate ener- 
gies (123, 148, 365, 402, 501, 697-702 and 943 cm"!) and 
five modes at high energies. We attribute these lorentzian 
lines to phonons. The factor group analysis yields 12 Ra- 
man active phonon modes according to the R — 3m space 
group: Tfiaman = 5 • Aig + 7 • Eg . In addition, we will 
later consider the A2g symmetry representations that can 
be induced by electronic resonances or DM interactionii. 
The corresponding tensors are: 

^a 0\ /c 0\ / e 

,Eg= -c d , A2g= -e 

^0 6/ \0 d Oj \ 

By rotating the samples within the ab plane (x — y 
plane) the Aig modes are expected to give a rotation 
invariant contribution in parallel (xx) polarization and 
disappear in crossed (xy) polarization of incident and 
scattered light. The Eg modes should show a four-fold 
symmetry in both polarizations. Due to the layered crys- 
tal structure and strongly anisotropic electronic polariz- 
ability perpendicular to the ab plane slight but inevitable 
misalignments of the crystals will show up as an addi- 
tional two-fold modulation of the scattering intensity. 

The 5 high energies phonons at about 3500 cm^^ (2 • 
Aig + 3 • Eg) are due to vibrations of hydrogen atoms. In 
a deuterated powder sample they shift to lower energies 
with a renormalization factor of 0.74. which is in good 
accordance to a simple estimate (Aw « ^/mii /"itt-d = 
0.71). This is depicted in Fig. 1 c) 
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We assign the modes at 365, 402 and 501 cm -^ to 



Eg symmetry, as their intensity clearly displays a four- 
fold rotational symmetry (see Fig. 1 b). The phonon 
modes at 123, 148 and 702 (697) cm^^ are distinctly 
different as they show a two-fold rotational symmetry. 
We assign these modes to Aig symmetry and the two- 
fold modulation of intensity to a small deviation of the k 
vector from the c axis, i.e. a slight canting of the sample 
surface. 

In a previous investigation of a series of 
Zn^Cu4_a;(OH)6Cl2 powder samples we stated a 
linear dependence of the phonon mode at 702 (697) 
cm^^ on composition. This allows us to estimate x = 0.8 
for the natural and a; = 1.0 for the synthesized single 
crystal, respectively. The synthesized single crystal is 
therefore in the regime of Herbertsmithite stoichiometry 
Still, the natural single crystal is reasonably close to this 
state in the sense that thermodynamic and spectroscopic 
experiments did not observe evidence for long range 
ordering for x = 0.8^. Studying the two different 
samples should therefore allow determining the influence 
of spin defects and static doping on the spin dynamics. 

A less intense, asymmetric mode with a fano lineshape 
is observed at 230 cm~^. This mode shows an even 
smaller intensity in the powder samples. Its lineshape 
may originate from the coupling of the corresponding 
phonon to a continuum of states. As its energy range 
is within the energy range of spin-fluctuationSi^i and the 
compound is an insulator the corresponding continuum 
is attributed to magnetic fluctuations. The mode itself 
is probably induced by crystallographic disorder as the 
larger intensity modes already exhaust the total number 
of Raman active modes. 

At low frequency and for high temperatures there is a 
pronounced quasi elastic scattering contribution (QES) 
visible with a two-fold rotational anisotropy (Fig. 1 b) 
that is similar to the phonon modes at 123, 148 and 702 
(697) cm~^, identified as Aig phonons. Therefore we as- 
sign the continuum to the same Aig component. At low 
temperatures (below T = 50 K) the quasi elastic scatter- 
ing is completely depressed as demonstrated in Fig. 2 for 
the case of the natural sample. 

An analysis of the QES (black dots) at 295 K is given 
in Figure 3 a) together with a Lorentzian (solid red line) 
and a Gaussian (dashed green line) fit. The data is in 
better agreement with the Lorentzian line shape. Figure 
3 b) shows the temperature development of the QES in- 
tensity in the synthesized (black dots) and the natural 
crystal (blue triangles). The decrease in intensity with 
decreasing temperature can be fitted with power laws 
(red and green curves). The observed exponent varies 
from I{uj ^ 0,r) - T3/2 to /(w ^ 0,r) - T^ for the 
synthesized and the natural sample, respectively. 

With decreasing temperatures we observe a moderate 
reduction of linewidth of all phonon modes, see Fig. 2 a. 
We attributed this to the decrease of anharmonic phonon 
scattering processes or other fiuctuations. The integrated 
intensity of the asymmetric Fano line shape at 230 cm~^ 
is increasing strongly; the intensity at 10 K is about 5x 
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FIG. 2. a) Temperature dependence of Raman spectra of 
the natural sample in the frequency regime 30 to 725 cm~^ 
and in parallel xx polarization. Spectra are shifted for clar- 
ity, b) and c) show the Bose corrected spectra for the nat- 
ural crystal and the synthesized sample, respectively, with 
phonon modes subtracted (dots) together with a fit to the 
background. The dashed black lines in b) and c) correspond 
to the fitted background in crossed xy polarization at 5 K. 
The spectra at T = 295 K in b) and c) are shifted in intensity 
for clarity. 



stronger compared to 120 K. In the synthesized sample 
this increase is only half that strong. At the same time, 
the linewidth of this mode is decreasing steadily before 
saturating for temperatures below approximately 70 K. 
With decreasing temperatures also the Aig phonons at 
123 and 148 cm~^ gain a Fano line shape. 

At low temperatures a broad, finite energy maximum is 
observed with increasing scattering intensity. This con- 
tinuum is observed for both the natural as well as the 
synthesized single crystal, as shown in Fig. 2 b) and 
c). For the synthesized sample the maximum position 
is higher {Emax = 265 cm~^ « 2.1J) compared to the 
natural one (215 cm^^ « 1.7J). In addition, the former 
sample also shows a more pronounced decrease of scat- 
tering intensity towards higher energies. We attribute 
the broad continuum in ZnCu3(OH)6Cl2 to scattering on 
magnetic correlations as it resembles two-magnon scat- 
tering of a quantum antiferromagnet in the paramagnetic 
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Further below we will discuss and compare its 
properties in detail with existing theories. The dashed 
lines in Fig. 2 b) and c) show the corresponding re- 
sults in crossed polarization. There are similar maxima 
and frequency shifts comparing the two sample qualities. 
However the onset of the scattering is shifted to higher 
energy for the natural sample. 

We observe a larger spectral weight of the continuum 
in the synthesized sample, especially in xx polarization. 
For both samples there is less spectral weight at low 
energies in xy polarization. This suggests the possibil- 
ity of decomposing the spectral weight into a low and a 
high energy part, to which the allowed symmetry chan- 
nels Aig, A2g and Eg contribute to varying extend. To 
estimate their magnitudes we investigate the rotational 
anisotropics of these signals in Fig. 5. It is important to 
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FIG. 3. a) Fit of Lorentzian (solid red line) vs. Gaussian 
(dashed green line) line shapes to the data (black dots) of 
the natural sample at low frequency and T — 295 K. The 
phonon modes at 123 cm~^ and 148 cni^^ are subtracted, 
b) Temperature dependence of the intensity of the QES in 
the synthesized crystal (black circles) and the natural sample 
(blue triangles). The dotted (dashed) line is a fit to the data 
using a power law according to / ~ y3/2 ^j ^ ys-j j-^j. ^j^g 
synthesized (natural) sample as described in the text. 



note that we observe no further signal that could be in- 
terpreted as of magnetic origin. In particular, there are 
no sharp modes at low energies. This statement is as- 
sured down to energy scales of approximately 5-10 cm~^ 
as no upturn is visible at the edge of the experimentally 
accessible energy window. 

Figure 4 a) zooms into the low frequency regime (up 
to 100 cm~^) at low temperature (T = 5 K). There is a 
linear frequency dependence towards the maximum and 
we can use its slope as a measurand of the integrated 
intensity, less sensitive to fluorescence backgrounds. In 
Fig. 4 b) the slope of the low energy scattering is plotted 
as function of temperature. It increases approximately 
linearly with decreasing temperature (T < 50 K). Both 
samples show this linear temperature dependence. How- 
ever at higher temperatures the synthesized sample has 
again the more pronounced T dependence as shown in 
Fig. 4 c). Such a linear increase of scattering intensity 
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FIG. 4. a) Linear fit to the low frequency scattering in xx 
polarization of the natural crystal at T = 5 K. b) Temperature 
development of its slope including a linear fit. c) Temperature 
development of the intensity of the continuum in both the 
natural and synthesized sample with guides to the eye. 
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FIG. 5. Angular dependence at T = 5 K of a) the low energy 
slope in the natural crystal and b) its continuum intensity. 
The data follow the rotational anisotropy expected for Eg and 
Aig symmetry components. Note that measurements were 
performed only in the range from 0° to 180°. 



is anomalous and not compatible with the expected Bose 
factor. However, it is frequently observed for magnetic 
Raman scattering in low dimensional spin systems that 
are close to quantum critical. 

Finally, in Figure 5 we probe the rotational anisotropy 
of the low energy slope in xx and xy polarization (shown 
in a) and the rotational anisotropy of the intensity of 
the broad continuum (shown in b) in the natural sample 
at r = 5 K. The low energy slope in a) shows a four 
fold Eg-like symmetry in parallel polarization, while in 
crossed polarization the rotational anisotropy is rather 
two- fold or of mixed character. The intensity of the broad 
continuum in b) shows slight anisotropics in both xx and 
xy polarization. 



IV. DISCUSSION 

In the following we will first discuss the QES observed 
at high temperatures followed by the implications of the 
low temperature data at intermediate and high energies. 

As mentioned above, the most characteristic feature 
of the QES is its strong depression of intensity with de- 
creasing temperatures. Such a scattering is frequently 
observed in compounds that realize low dimensional spin 
systems, as spin chains^ or frustrated dimer systems^!. 
It is related to energy density fluctuations that couple via 
spin-phonon coupling to the lattice and can be mapped 
on the specific heat Cm of the spin system, I{uj -^ Q,T) ^ 
Cm T^ ■ Its lineshape and symmetry should be Lorentzian 
and of Aig symmetry, respectively^^. The fit to the data 
given in Fig. 3 follows all these predictions very well. 
Again, the rotational invariance is distorted with two- 
fold symmetry but within experimental error identical to 
the one of the phonon modes identified as belonging to 
the Aig symmetry representation. 

From the temperature dependence of the quasi elastic 
scattering intensity /(w — >■ 0,r) we can derive informa- 
tion on the evolution of low energy excitations. In, e.g., 
the Shastry-Sutherland compound SrCu2(B03)2 the ex- 
ponential drop of the respective intensity has been used 
to derive a spin gaf>2i. The power law-like tempera- 
ture dependence in ZnCu3(OH)6Cl2 differs clearly from 
that and implies a spin liquid with more gradually evolv- 
ing correlations^^. We attribute the softer power law 
I{ijj — >■ 0, T) ~ T^/'^ in the synthesized single crystal to 
the reduced number of defects leading to an even more 
gradual evolution of spin correlations. 

We will now discuss the three most recent theoreti- 
cal investigations (using different approaches and tech- 
niques) of finite energy, magnetic Raman scattering in a 
kagome system and their implications. As was suggested 
by Cepas et alJ^, the formation of a valence bond crys- 
tal leads to a weakly broken translation symmetry of the 
spin system. Thereby the magnetic Raman response of 
the kagome plane should depend in a decisive way on the 
polarization of the incoming / scattered light, i.e. there 
should be characteristic oscillations of magnetic modes in 
the energy range of J when the polarizations are rotated 
within the plane. In contrast, the "true" spin liquid ex- 
hibiting no broken symmetry should show a polarization 
independent magnetic Raman response. Since experi- 
mental investigations using different techniques have up 
to now not observed a spin gap, the Raman experiment 
is highlighted to be potentially sensitive to even weakly 
broken symmetries in kagome plane. It should be noted 
that in the second order perturbational treatment used 
by Cepas et al. magnetic scattering is only observed in 
Eg symmetry, as the Ai^ component commutes with the 
used exchange scattering Hamiltonian. 

Recently, calculations based on a modified Shastry- 
Shrainian model describing a C/(l) Dirac spin liquid have 
been performed by Ko, et alJi. In this approximation 
also higher order terms, induced e.g. by electronic reso- 



nances of the phonon energy with correlated states of the 
system, are considered. Here, the included chiral s x s 
component leads to additional scattering contributions in 
Aig, A2C, and Eg symmetry components. From the Ra- 
man tensors given above we notice that the now allowed 
A2g component is only observed in xy polarization as a 
rotational invariant contribution. It does not contribute 
to XX polarization at all. The spectral distribution of 
the Raman scattering cross section derived by Ko, et al. 
shows a broad continuum that extends up to about 600 
cm^i (« 5- J). 

Finally, Lauchli et alJ^ calculated dynamical singlet 
fluctuations on the kagome lattice using large-scale exact 
diagonalizations on a 36 site lattice. The results show 
a broad continuum extending over a range of 2 — 3J. 
There is also a pronounced intensity shooting up at lowest 
energies, ut/J < 0.2 w 25 cm~^, caused by a high density 
of low energy singlet and triplet excitations. 

We can decompose the experimentally observed scat- 
tering continuum into two contributions both attributed 
to specific spin liquid excitations. The rotation depen- 
dence of the low energy slope (Fig. 5 a) shows in xx 
polarization a clear Eg-like contribution, while in the xy 
configuration it is oval shaped, most probably a mix of 
both Eg and A2g symmetry. The dominance of the Eg 
channel at frequencies below 100 cm~^ is in contrast to 
theoryil. The rotation dependence of the overall con- 
tinuum's intensity (Fig. 5 b) shows a weak asymmetry 
in both XX and xy configurations for the natural and 
the synthesized sample. This anisotropy has no domi- 
nant four-fold contribution. However, a superposition of 
a four-fold and a two-fold contribution leads to a good 
description of the scattering intensity. We conclude that 
the Aig (for xx) and the A2g (for xy) symmetry channels 
are major contributions for the continuum at an energy 
of « 2J. These components are proposed to be due to 
higher order perturbation theorj»ii. Comparing our ex- 
perimental data in crossed polarization to the data in 
parallel polarization, no clear separation of the respec- 
tive spectral weight is evident. 

This behavior is in clear contrast to the phonon Ra- 
man scattering, where modes of Aig and Eg show a pro- 
nounced selection rule and different in-plane anisotropy. 
The observed Fano-like line shape of phonons is also only 
observed for modes with Aig symmetry. On the other 
hand phonon and magnetic scattering differ within two 
respects that in the theoretical literature on the Her- 
bertsmithite are jointly discussed^. One is related to 
defects, which for the local spin excitations play a dif- 
ferent role than for g w optical phonons with large 
coherence length. The second aspect is DM interaction, 
which mixes singlet and triplet excitations. DM interac- 
tion and structural defects as intersite mixing will allow, 
further enhance and possibly mix Aig and A2g contribu- 
tions, while spin defects in the plane reduce the scattering 
intensity. Therefore, a clear assignment of contributions 
from the Aig, A2g and Eg channels to the background 
is not straight forward. Still, the energy range and gen- 



eral shape of the observed background are in very good 
agreement with the above mentioned calculations of Ko, 
et al.il. 

The effects of impurities in our Herbertsmithite sam- 
ples should be considered, and these effects should more 
pronounced for the natural sample (x = 0.8). While 
keeping the kagome plane intact, Cu atoms in the in- 
termediate layer may induce lattice deformations due 
to Jahn- Teller distortion, as well as additional magnetic 
coupling to the kagome Cu atoms. The effect of impu- 
rities in quantum spin systems with dedicated low en- 
ergy excitations can be exemplified using the frustrated 
dimer system SrCu2 (603)2, where f-2% of Zn doping 
on Cu broadens the low energy excitations and consider- 
ably redistributes their spectral weighli^^. This may be 
the reason that in contrast to theory, our experimental 
data do not show sharp magnetic modes or a pronounced 
low energy maximum. In NMR experiments on Herbert- 
smithite it has been shown that the local susceptibility 
in the proximity to a defect is different and anomalous 
compared to the regular behavior—. With respect to Ra- 
man scattering it can be assumed that the excitations in 
the continuum that consist mainly of two-pair statesii) 
should be reduced in intensity, comparing the natural 
and the synthesized sample. Such a reduction of contin- 
uum scattering intensity is indeed observed, as well as a 
softening of about 20%. 



V. CONCLUSION 

Our Raman data give further evidence that the alge- 
braic spin liquid is the ground state of the Heisenberg 



s = 1/2 kagome lattice antiferromagnet. This is based 
on the observation of a power law depletion of thermally 
induced fluctuations of the magnetic energy density with 
decreasing temperature and scattering continua in the 
low- and mid energy range. We give estimates of the con- 
tributions to these scattering continua to different sym- 
metry components and conclude that except for the low 
energy continuum in xx configuration, no strong Eg con- 
tribution has been observed. 

Our data leads to the conclusion that the kagome plane 
has short range and continuous correlations that form a 
maximum (at w 2 J) with no clear cutoff but a moderate 
suppression at higher energy scales (at w 6J). The over- 
all energy range and shape of the continua agrees well 
with existing calculations for a spin liquid state. 

Presently the role of impurities in the single crystal, i.e. 
antisite disorder and local, defect-like interlayer interac- 
tions, is difficult to evaluate. Comparing the cleaner syn- 
thesized with the natural, natural samples a weak sup- 
pression of the continuum intensity with a frequency shift 
is evident. As no sharp, pronounced magnetic peaks are 
observed their relation to disorder remains an open issue. 
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